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ABSTRACT 

A method is presented to model atomic oxygen erosion of protected polymers 
in low Earth orbit (LEO). Undercutting of protected polymers by atomic oxygen occurs 
in LEO due to the presence of scratch, crack or pin-window defects in the protective 
coatings. As a means of providing a better understanding of undercutting processes, a 
fast method of modeling atomic -oxygen undercutting of protected polymers has been 
developed. Current simulation methods often rely on computationally expensive ray- 
tracing procedures to track the surface-to-surface movement of individual “atoms”. 
The method introduced in this paper replaces slow individual particle approaches by 
substituting a model that utilizes both a geometric configuration-factor technique, 
which governs the diffuse transport of atoms between surfaces, and an efficient 
telescoping series algorithm, which rapidly integrates the cumulative effects stemming 
from the numerous atomic oxygen events occurring at the surfaces of an undercut 
cavity. This new method facilitates the systematic study of three-dimensional 
undercutting by allowing rapid simulations to be made over a wide range of erosion 
parameters. 


1. INTRODUCTION 

Undercutting of protected polymers by atomic oxygen occurs in LEO due to 
the presence of scratch, crack or pin-window defects in the protective coatings. 
Methods to accurately model characteristics of undercutting beneath such defects exist, 
but due to the extremely complicated nature of the problem fall short of providing a 
complete representation of the undercutting process. For example, Monte-Carlo 
methods [1-2], when used in conjunction with ray-tracing procedures, produce 
remarkably faithful representations, yet are computationally exhaustive for small cell 
sizes, and thus often restricted to two spatial dimensions. Another deficiency for any 
method is the lack of knowledge as to what effective values should apply for atomic 
oxygen reaction and recombination rates at polymer and protective coating surfaces. 
Such lack of knowledge complicates correlation between in-space and ground based 
results. 
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Proper correlation between ground-based systems and flight results are 
required for meaningful durability studies. A means of establishing such a correlation 
is the use of computer codes to simulate ground-based and in-space environments. To 
advance modeling of three-dimensional undercutting erosion yet minimize 
computational requirements, a procedure has been developed to provide fast and 
accurate simulations for a wide variety of conditions. This method shortcuts many 
computational bottlenecks and avoids lengthy ray tracing techniques by applying 
geometric configuration factors to determine the exchange of atoms between cavity 
surfaces, which in turn are approximated by simple geometric shapes. 

The primary goal was to construct a procedure that provides reasonably 
accurate predictions and serves as a guide for the design of more complex codes. 
Simple geometries that are easily parameterized were selected. To limit the scope, 
individual geometries were chosen that exhibit qualities peculiar to directed beam 
undercutting resulting from a fixed ram direction. Therefore the geometric models used 
in this paper do not apply to simulation of arbitrary three-dimensional cavities, but are 
restricted to modeling classes of cavities that exhibit symmetry about an axis or a plane. 


2. GEOMETRIC MODEL 

As mentioned above, the code models three-dimensional undercut cavities by 
using a few prescribed geometries to serve as representative shapes for some of the 
more simple cavities observed from in-space undercutting below pin-window defects. 
In the model, undercutting occurs directly below a defect present in the upper protective 
coating of the film. With a double-coated film, a lower protective coating also exists. 
The protective coating is assumed to be infinitely thin. It is convenient and reasonable 
to use circular disks to represent the pin-window defects. Consistent with this simple 
defect shape, all model-geometry cavity cross sections in planes parallel to the polymer 
film are circular. In the present code, four specific geometric classifications exist. 
These four geometries are: (1) normal cylinder, (2) normal-truncated cone, (3) oblique- 
cylinder, and (4) oblique-truncated cone. Each of these individual models was chosen 
to capture a specific set of features characteristic of the posthole-type erosion pattern 
arising from undercutting due to directed beam atomic oxygen. The degree of 
obliqueness of model (3) or model (4) is determined by the angle that the ram direction 
deviates from the film’s surface-normal direction. The conical geometries allow for 
modeling of converging and diverging undercut cavities. 

An equation for volume that is valid for the above geometric models is given 
by the following equation: 


V =J h ( r i 2 + r i r 2+ r 2 2 )' (!) 

where r ; and r 2 are the upper and lower circular disk radii (see Figure 1), which are 
equal in the case of a cylindrical model. In order to provide a reasonable 
characterization, the radial and vertical dimensions of the model cavity must evolve 
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synchronously with cavity volume in a suitable fashion. To automate the process and 
provide a degree of control in specifying model geometries, the radial coordinate, r u is 
chosen to be a simple function of the cavity height, h, for h < /; k , where h k is the 
polymer-film thickness. This function is given as follows: 

t\ = r def + ah , h < h k . (2) 

For cylindrical geometries the radius is initially equal to the defect radius, r de f, and 
increases linearly with cavity height for h < h k . In the case of conical geometries an 
additional parameter, b, is introduced to relate the radii of the circular disks. Equation 
(2) is again used to specify r k , and then r 2 is put in terms of r k as given by the following 
relation: 


r 2 =r l +bh, h<h k . (3) 

As the value of the parameter b is negative, zero, or positive, the conical section is 
converging, cylindrical, or diverging, respectively. Since the volumes of these models 
are known functions of the radius and height, the cavity height can be solved in terms 
of the volume. The cavity’s height, h, cannot exceed the polymer-film thickness, h k . 
The protective-coating thickness (typically 0.13 micron) is ignored. In this manner the 
cavity geometry is specified once the cavity volume is known. The parameters a and b 
control the cavity aspect ratio and cavity wall slope, respectively. As mentioned above, 
the degree of variation in the ram direction from normal determines the obliqueness of 
the model. 


3. PHYSICAL MODEL 

The present model assumes that non-reacting atomic-oxygen atoms reflect off 
surfaces diffusely, and consequently their movements between surfaces can be treated 
mathematically like diffuse radiation exchange between surfaces. The fraction of 
uniform diffuse radiation leaving one surface that reaches another is equivalent to the 
configuration factor between two surfaces, because it depends solely on the geometric 
orientation of the surfaces with respect to one another [3]. The geometric dependence 
of configuration factors can be used to derive algebraic relationships between factors. 
One pair of configuration factors exists for each pair of finite surfaces. A reciprocity 
relationship exists for the factors between two finite surfaces. It can be expressed by 
the following formula: 


A F : 


I '~J 


A /h’ 


(4) 


where 
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(5) 


^4 


M f 

d J A J/ 


cos 0, cos 0, 

-z — - <r/A . <rZA ; 

'a, ;r 5 2 7 ' 


is the configuration factor from area Aj to area A,. The angles Q. and Q . are the angles 
between the line segment S, connecting the differential areas dA, and dAj , and surface 
normals n, and /i ; , respectively. To complete the general relations among configuration 
factors needed here, a conservation relation stating that the sum of fractions of emission 
leaving a surface and arriving at other surfaces of an enclosure (including the emitting 
surface) must sum to unity is expressed by equation (6). With the emitting surface 
denoted by the subscript index j and the N receiving surfaces of the enclosure identified 
by the running subscript index i, then the sum of configuration factors can be written 
compactly as: 


1 ^ = 1 . «) 

i 


Given the above set of relationships, a complete set of configuration factors 
can be derived for an enclosed system starting from a relatively small initial subset, 
which for the geometries used here can be obtained from the existing literature on 
configuration factors [3], An oblique truncated cone (converging downward from the 
defect surface) is shown in Figure 1. This geometry is one of the simple cavity 
geometries used here for thin film modeling. The upper surface of the enclosure is 
comprised of a central disk that represents the defect, and a circular annulus that 
represents that portion of the coating from which the substrate has been removed by 
undercutting. These two areas are denoted as area 1 and area 2, respectively. 
Connecting the upper and lower disks is the lateral wall area denoted as area 3. The 
lower disk is denoted as area 4. Assuming F 4 _ i and F 4 . ( i+ 2 ) are either known by 
catalogued formula or calculated directly using equation (5), where the summed 
subscripts in the last factor indicate that the corresponding surface is a composite 
surface, then the complete set of factors for the enclosure is found by simple algebra. 
A list of this set of configuration factors and their origin are given in Table 1. 


4. COMPUTATIONAL MODEL 

The first step in arriving at a computational model is to provide an algorithm 
to determine the fraction of atoms arriving at a given surface after surviving an arbitrary 
number of bounces within the enclosure. Let f° be the fraction of atoms arriving 
directly through the defect that reach the i surface without reflecting off a surface, and 
let q! be the probability that an atom arriving at the i surface after n previous cavity 
bounces will not react, recombine, or exit during the next surface collision, but reflect 
again. The fraction of atoms surviving one reflection is q'f-' . Thus, by summing over 
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all reflecting surfaces j=l, 2,3... N and employing the configuration factor, the fraction 
of atoms reaching the i th surface after one bounce is given by the following sum: 

fi' 4° f! ■ < 7 > 

j 

Although the configuration factors must be recalculated as the cavity grows, they are 
held stationary over a given time-step interval. By telescoping from the first bounce, it 
can be shown that the fraction, f ." , of initial atoms entering the cavity during a given 

time-step interval that reach the i surface after n bounces is given by the following 
summation over the N surfaces: 

= " = 1.2,3,... (8) 

j 

The fraction of atoms recombining at a polymer surface, f p , is found by 
summing (over the total number of impacts in a time step) the product of the fraction of 
atoms arriving at the surface times the probability P" of reaction, which in general 
changes between surface impacts. On average, energy is lost following a collision, and 
this loss is generally reflected in the code by selecting a lower reaction probability for 
the next impact. For convenience, consider the set of all polymer surfaces as being one 
surface denoted by the subscript p. Then the net fraction of atomic-oxygen atoms 
reacting at the polymer to cause erosion during a time step consisting of NB bounces is 
given by the following expression: 


NB N 

/r=s^'iA,c‘/r'. (d 

H=1 j 

The total number of atoms reacting is given by the product of the net fraction 
of atoms reacting and the total atoms entered, which itself is given by multiplying the 
atom fluence, A F , during a time-step interval by the defect area, A,/,./. Continuing in 
this fashion, the product of the total number of atoms reacting and the erosion yield, E, 
(volume/atom) equals the amount of volume erosion, AV , obtained during a time step, 
which is given by the following: 


AV = /JA FA def E. 


( 10 ) 


The final undercut volume, V , is obtained by computing a series of AV volumes. 

In preparation for discussing results, it is convenient to define the volume 
erosion gain, G, as the ratio of the volume of undercutting erosion below a defect site to 
the volume of erosion obtained on the surface of an unprotected “smooth-witness” 
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sample, which had been exposed to the same number of atoms that entered the defect. 

the above gain can be determined by dividing the 
covering N s time steps by the initial reaction 

X rNB 

n — . ( 11 ) 

P> S 

The “instantaneous” gain, G ' , corresponding to the gain for a single time-step interval 
is given by 

G' = f™IP% ( 12 ) 

and provides a normalized erosion rate useful for comparison purposes. When 
G' equals one, the rate of undercutting erosion equals the rate of erosion of an 
unprotected smooth sample of area A de f exposed to the same atom flux. 


In terms of simulation parameters, 
average of the incremental f 
probability, P° , and is given by 


5. RESULTS AND DISCUSSION 

Values of the erosion parameters representative of the LEO environment were 
used in simulations. These nominal values along with other code parameters are listed 
in Table 2. Any variation from these nominal values is noted in the text as the 
case occurs. Similar to ground-based experiments, simulated fluence is an effective 
fluence, equivalent to the atoms per unit area in LEO that would produce the same 
volume erosion. The in-space erosion yield E for Kapton® representative of LEO is 
3 x 10' 24 cm 3 /atom based on atomic oxygen energy of 4.5 eV. It is emphasized that the 
tabulated entry for atomic-oxygen fluence is scaled by polymer-film thickness. Thus, 
normalized fluence-per-time-step-interval is listed as 2X 10 22 atoms/cm 2 per cm of film 
thickness, giving a non-scaled cumulative fluence over 100 time steps, for example, of 
5 X 10 21 atoms/cm 2 for a 0.0025 cm thickness film. In this fashion, having normalized 
the erosion rate by film thickness, the results in this paper hold for any film thickness. 

To give an example of effects obtained by changing cavity shape, results are 
given comparing data for the normal-cylinder cavity model with data for the converging 
and diverging normal-cone cavity models. For the conical geometries, the final wall 
angle converges or diverges (downward from the defect) by 7.1 degrees from axial. For 
these three geometric models, the variation of normalized erosion rate G' versus 
fractional fluence is presented in Figure 2. It is readily seen in Figure 2 that whereas 
each curve is distinguished by the location of a steep dip, the characteristic shape of the 
curves remains the same for the different geometries. The dip in a curve coincides with 
the cavity reaching the lower protective coating. The reason that the dip occurs at a 
lower fluence for a converging cavity is that its volume grows more slowly with cavity 
height than the other cases, causing it to reach the lower coating first. For the same 
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reason, the cavity of the cylindrical model reaches the lower coating before the cavity 
of the diverging-cone model. 

Two of the most important parameters influencing undercutting erosion are the 
defect size r de f and reaction probability P' 1 . To identify any gross effects observed with 
changes in these two parameters, the variation of final erosion gain, G, as a function of 
r def is shown in Figure 3 for various values of P° using the normal-cylinder geometry. 
For this data, r def ranges from 0.2 h k to 0.001 h k . For a given P°, G increases rapidly 
with decreasing r def until, at a value of r def approximately equal to 0. 1 h k , little increase 
in G is observed. Gains for the largest defect size, r de f = 0.2 h are close to one, while 
gains for the smallest defect size, r de f = 0.001 h k , are about 75% greater because of 
enhanced trapping of the atomic oxygen. It is seen that variations of plus or minus 20% 
from the nominal reaction probability ( P° = 0.11) produce respective changes in G of 
approximately -2.1% and 2.2% for the smallest defects and -0.6% and 1.5% for the 
largest defects. 

As a final example, a comparison of erosion rates is presented for single- 
coated and double-coated Kapton® films. To simulate the effect of having no lower 
protective coating, the recombination probability is set equal to one at the lower film 
surface. In Figure 4 curves of the erosion rate, G' , as a function of fractional fluence 
are shown for single -coated and double-coated Kapton® film using the normal-cylinder 
model. As expected, the curves are identical until the bottom surface is reached. In the 
double-coated polymer case, the presence of the lower coating results in additional 
erosion, beyond that obtained in the single-coated case, due to reflected atoms. 


6. SUMMARY 

A new method to simulate three-dimensional undercutting below pin-window 
defects in protected polymer thin films has been constructed. This method is very fast. 
To minimize computational requirements, it uses a geometric configuration-factor 
technique to govern the exchange of atoms between surfaces. The model assumes that 
non-reacting, non-recombining atoms reflect diffusely off surfaces. The basic 
equations to derive the view factors are given. In addition, a table listing the sequence 
of steps to derive a full set of configuration factors is provided for a representative case. 
A set of simple geometric shapes is used to model undercut cavities. A given geometric 
shape is selected based on its suitability to faithfully model particular qualities. These 
qualities belong to a small set of cavity characteristics, which include the cavity aspect 
ratio, cavity obliqueness to the film surface due to ram angle, and cavity-wall 
divergence. The equations relating the amount of volume erosion and the rate of 
volume erosion to the atomic-oxygen fluence are presented. Using a suitable set of 
simple geometries as models, rapid simulations can be obtained over a wide range of 
parameter space. Examples are given that illustrate typical results that may be obtained 
using the associated computer code. 
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7. TABLES 


Table 1. List of Configuration Factors Associated with the Geometry Illustrated in Figure 1 and Their Origin 


Configuration Factor 

Origin of Configuration Factor 

F 4 -i 

Known 

E 4-(\+2) 

Known 

Fa-2 = E a-(\+2)~E A-\ 

From definition of F4-(i+2> = F4-1 + Fa -2 

7^ 

11 

11 

S 

L 

II 

O 

Planar Areas, Equation (5) 

O 

II 

£ 

Area 1 and area 2 in same plane, equation (5) 

E\-a — Aa Ea-\/A\ 

Reciprocity, equation (4) 

F 2 -A = Aa F4-2/A2 

Reciprocity, equation (4) 

F 1-3 — 1- F 1-4 

Conservation of atoms from area 1 , equation (6) 

F2-3 = 1 - F 2 -a, 

Conservation of atoms from area 2, equation (6) 

^4-3 = 1- F4-I- F4.2 

Conservation of atoms from area 4, equation (6) 

F3-1 = A\ F1-3/A3 

Reciprocity, equation (4) 

F3-2 = A2 F2-3/A3 

Reciprocity, equation (4) 

F3-4 = A4 F4.3/A3 

Reciprocity, equation (4) 

t 

tC 

11 

£ 

Conservation of atoms from area 3, equation (6) 


Table 2. Nominal Computational Model Parameters for LEO Atomic Oxygen Interaction with Kapton ( 


Atomic oxygen initial impact at polymer reaction probability 

0.11 

Atomic oxygen asymptotic impact at polymer reaction probability 

0.001 

Atomic oxygen impact at coating recombination probability 

0.13 

Atomic oxygen impact at polymer recombination probability 

0.24 

Fractional decrease in polymer reaction probability per impact 

0.368 

Number of atom bounces per time step 

500 

Atomic oxygen fluence (atoms/cm 2 ) per cm of film thickness per time step 

2X 10~ 

Number of erosion time steps 

100 

Atomic oxygen fluence (atoms/cm 2 ) for 0.0025 cm thickness film 

5 X 10 21 

Erosion yield (cm 3 /AO-atom) 

3 X 10' 24 

Value of parameter a in equation [2] 

0.0759 

Values of parameter b for h = h^ in equation [3] for conical geometry 

1+ 

A 

Protective coating defect size 

0.1 /lie 
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Figure 3. 
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